Heralded Entanglement between Atomic Ensembles: 
Preparation, Decoherence, and Scaling 
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Heralded entanglement between collective excitations in two atomic ensembles is probabilistically 
generated, stored, and converted to single photon fields. By way of the concurrence, quantitative 
characterizations are reported for the scaling behavior of entanglement with excitation probability 
and for the temporal dynamics of various correlations resulting in the decay of entanglement. A 
lower bound of the concurrence for the collective atomic state of 0.9 =b 0.3 is inferred. The decay of 
entanglement as a function of storage time is also observed, and related to the local dynamics. 
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Beyond a fundamental significance, quantum control 
of entanglement between material systems is an essential 
capability for quantum networks and scalable quantum 
communication architectures [3,01. In recent years, sig- 
nificant advances have been achieved in the control of 
the quantum states of atomic systems, including entan- 
glement of trapped ions [3|, |4( and between macroscopic 
spins By following the seminal paper of Duan, Lukin, 
Cirac and Zoller (DLCZ) [6], entanglement between sin- 
gle collective excitations stored in two remote atomic en- 
sembles has also been demonstrated 0- In the DLCZ 
protocol, entanglement is created in a probabilistic but 
heralded way from quantum interference in the measure- 
ment process d, 0, [To| . The detection of a photon from 
one or the other atomic ensemble in an indistinguishable 
fashion results in an entangled state with one collective 
spin excitation shared coherently between the ensembles. 
Such entanglement has been critical for the initial imple- 
mentation of functional quantum nodes for entanglement 
distribution for the investigation of entanglement 

swapping [12[ and for light-matter teleportation [l3[. 

Because of the relevance to quantum networking tasks, 
it is important to obtain detailed characterizations of the 
physical processes related to the creation, storage, and 
utilization of heralded entanglement. Towards this end, 
significant advances have been demonstrated in the gen- 



of collective excitation 
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eration of photon-pairs 14 . 1 1511 and the efficient retrieval 



Moreover, decoherence 



processes for a single atomic ensemble in the regime of 
collective excitation have been investigated theoretically 
[lit and a direct measurement of decoherence for one 
stored component of a Bell state recently performed [19[ . 
However, to date no direct study has been reported for 
the decoherence of an entangled system involving two 
distinct atomic ensembles, which is a critical aspect for 
the implementation of elaborate protocols [13, E|, |22| ]. 
The decoherence of entanglement between ensembles has 
been shown in recent setups, through the decay of the 
violation of a Bell inequality I 111 and the decay of the 
fidelity of a teleported state [l3[. However, a quantita- 
tive analysis was not provided since these setups involved 



many others parameters, such as phase stability over long 
distances. 

In this Letter, we report measurements that provide 
a detailed and quantitative characterization of entangle- 
ment between collective atomic excitations. Specifically, 
we determine the concurrence C [23[ as a function of 
the normalized degree of correlation gi2 [l6[ fc> r the en- 
sembles, including the threshold for entanglement 
(C > 0). We also map the decay of the concurrence C(r) 
as a function of storage time for the entangled state, and 
interpret this decay by measuring the local decoherence 
on both ensembles taken independently. Compared to 
Ref. [7], these observations are made possible by a new 
system that requires no active phase stability and that 
implements conditional control for the generation, stor- 
age, and readout of entangled atomic states. 

Our experiment is illustrated in Fig. [H A single cloud 
of cesium atoms in a magneto-optical trap is used; two 
ensembles are defined by different optical paths 1 mm 
apart [ll|, Q . This separation is obtained by the use of 
birefringent crystals close to the cloud, which separate 
orthogonal polarizations [25[. At 40 Hz, the trap mag- 
netic field is switched off for 7 ms. After waiting 3 ms 
for the magnetic field to decay, the two samples are si- 
multaneously illuminated with 30 ns-long and 10 MHz 
red-detuned write pulses, at a rate of 1.7 MHz. Given 
the duty cycle of the experiment, the effective rate is 180 
kHz. Spontaneous Raman scattered fields induced by the 
write beams are collected into single- mode fibers, defin- 
ing for each ensemble optical modes that we designate 
as fields 1u,d with 50 jim waist and a 3° angle relative 
to the direction of the write beams [HI Gil- Fields 1u,d 
are frequency filtered to block spontaneous emission from 
atomic transitions \e) — > which do not herald the 
creation of a collective excitation. After this stage, and 
before detection, fields 1u,d are brought to interfere on 
a polarizing beam-splitter. A detection event at £>i a ,i6 
that arises indistinguishably from either of the fields 1u,d 
projects the atomic ensembles into an entangled state 
where, in the ideal case, one collective excitation is co- 
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FIG. 1: (a) A weak write pulse is split into two paths sepa- 
rated by 1 mm and excites simultaneously two atomic sam- 
ples, U, D. The resulting fields 1u,d are combined at the 
polarizing beamsplitter (PBS) and sent to the single-photon 
detectors ^i a ,i6- A detection event at D\ a or D\b heralds the 
creation of entanglement, (b) After a storage time t, entan- 
glement is verified by mapping the atomic state to propagat- 
ing fields 2u,d by way of read pulses. Tomography is then 
achieved in two steps, as described in the text. The atomic 
cloud is initially prepared in the ground state \g). {\g), |e)} 



denote the levels {\6S 1/2 ,F = 4), \6S 1/2 ,F = 3), |6R 



3/2 



F 



4)}. Note that the fields 1u,d and 2u,d are detected with a 
small angle relative to the classical beams, which is not rep- 
resented here for the sake of simplicity. 



herently shared between the J7, D ensembles 

In the ideal case of small excitation probability, the 
atom- field 1 joint state can be written for each ensemble: 



|*) = |o )|0i) + vx|ia)|ii> + o(x) 



(1) 



with \m) the state of the field 1 with n photons and \n a ) 
the state of the ensemble with n collective excitations. 
Upon a detection event at £>ia,i6, in the ideal case, the 
atomic state is projected into 



U,D 



V2 



{\Va)u\la)D±e ie \l a )u\Qa)D)+0(Jx), (2) 



where \O a )u,D, \la)u,D refers to the two ensembles Z7, D 
with 0,1 collective excitations, respectively [6]. The ± 
sign is set by the detector that records the heralding 
event. The overall phase is the sum of the phase differ- 
ence of the write beams at the U and D ensembles and the 
phase difference acquired by fields 1 in propagation from 
the ensembles to the beamsplitter. To achieve entangle- 
ment, this phase must be constant from trial to trial [261 ]. 
In order to meet this requirement, the initial demonstra- 
tion reported in [7[ employed auxiliary fields to achieve 
active stabilization for various phases for two ensembles 
located in distinct vacuum apparatuses. By contrast, in 
our current setup (Fig. QJi), is determined only by the 
differential phase for the two paths with orthogonal po- 
larizations defined by the birefringent crystals [25]; our 
small setup has sufficient passive stability without need 
of adjustment or compensation as the phase does not 
change by more than a few degrees over 24 hours. 



To verify operationally entanglement between the U, D 
ensembles, the respective atomic states are mapped into 
photonic states by applying simultaneously read pulses in 
the configuration introduced in Ref. [15[, as depicted in 
Figure [TJd. The delocalized atomic excitation is retrieved 
with high efficiency thanks to collective enhancement 

[l6[ and, in the ideal case, \^u,d) would be mapped 
directly to the photonic state of fields 2u,d with unity ef- 
ficiency and no additional components. Stability for the 
phase difference of the read beams and of fields 2u,d is 
also required in this process; it is again achieved by the 
passive stability of our current scheme [25]. Since entan- 



glement cannot be increased by local operations [27[ , the 
entanglement for the atomic state will always be greater 
than or equal to that measured for the light fields. 

A model-independent determination of entanglement 
based upon quantum tomography of the fields 2u,d has 
been developed in Ref. |7[. The model consists of re- 
constructing a density matrix, p2 Vl 2 D ^ obtained from the 
full density matrix by restriction to the subspace with no 
more than one photon per mode. We also assume that all 
off-diagonal elements between states with different num- 
bers of photons vanish. The model thus leads to a lower 
bound for entanglement. As detailed in Ref. Q, p2 v ,2 D 
can be written in the photon- number basis \n)\m) with 
{n, m} = {0, 1} as follows: 



P2u,2d = d 



/poo \ 

poi d 

d* p 10 

V Pll J 



(3) 



Here, pij is the probability to find i photons in mode 2jj 
and j in mode 2]j; d is the coherence term between the 
|1)|0) and |0)|1) states; and P = poo + Poi + Pio +Pn- 
From p2u,2 D , one can calculate the concurrence C, which 
is a convenient monotone measurement of entanglement 
ranging from for a separable state to 1 for a maximally 
entangled state [231 ] : 

C = max(0,C ) with C = ^(2\d\ - 2^/mpH). (4) 

In the regime of low excitation and detection probabilities 
in which the experiment is performed, the vacuum poo 
can be approximated by poo ~ 1 — p c , while the terms 
poi and p 10 are given by p 10 = p 01 = p c /2. p c is the 
conditional probability of detecting a photon in field 2 
from one ensemble following a detection event for field 1. 

Experimentally, we reconstruct p2u,2 D and then calcu- 
late C by using two configurations for the detection of 
fields 2 corresponding to two settings of the (X/2) v 
waveplate shown in Fig. [TJd. The diagonal elements of 
P2 Uy 2 D are determined from measurements of the photon 
statistics for the separated fields 2^/, 2^, i.e., by detecting 
independently each field. To access the coherence term 
<i, fields 2jj d are coherently superimposed and the count 
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FIG. 2: Concurrence C as a function of the normalized cross 
correlation function gi2, for the two possible heralding events 
(detection at D\ a or D\b). The solid line is obtained from 
Eq. (|5]) with the fitted overlap (see inset) and assuming an 
independently-measured retrieval efficiency at 13.5%. The 
dotted line corresponds to Co. Inset: Average visibility of the 
interference fringe between the two field- 2 modes. The solid 
line is a fit using the expression given by Eq. (0, with the 
overlap f fitted to 0.95 ± 0.01. 



rates from the resulting interference are recorded as a 
function of the relative phase between the 2u,d fields. It 
can be shown that d ~ V(pio+Poi)/2 ~ Vp c /2 [7J, where 
V is the visibility of the interference fringe. 

To investigate the scaling of entanglement with exci- 
tation probability we determine C for various values 
of x f° r fixed memory time r = 200 ns. As \ increases, 
higher order terms in the expansion of Eq. ([2|) cannot be 
neglected, precisely as in parametric down conversion. A 
convenient parameter to assess the excitation regime of 
each ensemble is the normalized intensity cross correla- 
tion function gyi between field 1 and field 2 [l6[, defined 
as #12 = P12KP1P2) with P12 the joint probability for de- 
tection events from field 1 and 2 in a given trial and pi 
the probability for unconditional detections in field i. In 
the ideal case, this function is related to the excitation 
probability x by gi2 = 1 + 1/x, where #12 > 2 defines 
the nonclassical border in the ideal case [14| and gi2 ^> 2 
being the single-excitation regime for the ensembles. 

Expressing the two-photon component for the two en- 
sembles as pn = xPc 00 Pel ^12 ? we rewrite the concur- 
rence as: 



C ~ max[0,p c (V - 2^(l-p c )/g 
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(5) 



a {D) 



where g±2 is for either ensemble alone, with g\ 2 } 
g\2 assumed. The visibility V can be expressed in terms 
of gi2 as the higher order terms act as a background noise. 
With (l/2)pi£>2 a good estimation for the background, 
the visibility can be written as [25[ 



P12 ~ P1P2 
P12 +P1P2 
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where £ is the overlap between fields 2u,d [2S|- In the 
limit of near zero excitation, as gi2 goes to infinity, the 
concurrence reaches its asymptotic value given by the 
retrieval efficiency ^p c (29[. 

Fig. [2] presents our measurements of the concurrence 
C as a function of #12. As the excitation probability is 
decreased, gi2 increases as does the entanglement. The 
threshold to achieve C > is found to be g$ — 7, corre- 
sponding to a probability p ~ 1.2 x 10 -2 per trial for the 
creation of the heralded entangled state and to a prepa- 
ration rate ~ 2 kHz. Note that C = (or C not greater 
than zero) does not imply that there is no entanglement, 
only that any possible entanglement is not detected by 
our protocol, which provides a lower bound for the en- 
tanglement. More importantly, in an infinite dimensional 
Hilbert space, entangled states are dense in the set of all 
states [30| , so that zero entanglement is not provable for 
an actual experiment by way of the concurrence. 

To confirm the model leading to Eq. (|6]), the inset gives 
the measured visibility V as a function of gw The solid 
line is a fit according to Eq. (|6|) with free parameter £, 
leading to an overlap £ = 0.95 ±0.01, in agreement with 
the value £ = 0.98 ± 0.03 obtained from an independent 
two-photon interference measurement. With the fitted 
value of £ and with the independently determined value 
of the conditional probability p c = 0.135 ± 0.005 from 
measurements performed on each ensemble separately, 
we compare our measurements of C with the prediction 
of Eq. [5] (solid line in Fig. [2j) and find good agreement. 

Table I provides the diagonal elements of the den- 
sity matrix p2 Vl 2 D and the concurrence for the case 
gi2 = 60 ± 4 corresponding to a probability to create 
entanglement p = 9 x 10 -4 per trial (160 Hz). A value 
C = 0.092 ± 0.002 is directly measured at detectors 
D2a,D2b without correction. By way of the indepen- 
dently determined propagation and detection efficiencies, 
we infer the density matrix for fields 2^/, 2^> at the 

output of the ensembles, from which we obtain a con- 
currence c% ut n Ut = 0.35 ± 0.1. This value exceeds the 
published state of the art by two orders of magnitude [7|] . 
This leap underlines the progress obtained in terms of 
suppression of the two-photon component and achievable 



TABLE I: Diagonal elements and concurrence of the density 
matrices for fields 2u,d, without and with correction for prop- 
agation losses and detection efficiencies. The last column pro- 
vides the estimated elements and concurrence for the atomic 
state by considering the readout efficiency 77. g±2 = 60 =b 4. 
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(6) 







-output 


PU,D 


Poo 


0.864 ±0.001 


0.54 ±0.08 


0±0.3 


Pio 


(6.47 ±0.02) x 10- 2 


(22 ±4) x 10-2 


0.5 ±0.15 


P01 


(7.07 ±0.02) x 10- 2 


(24 ±4) x 10-2 


0.5 ±0.15 


Pll 


(2.8 ±0.2) x 10- 4 


(3 ±2) x 10- 3 


0.015 ±0.025 


C 


0.092 ±0.002 


0.35 ±0.1 


0.9 ±0.3 
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0.05 



0.00 




decay. The temporal dynamics reveal a finite-time decay, 
with separability onset for storage time r ~ 20 jis. From 
a more general perspective, the inferred concurrence for 
the collective atomic state, Cu,d = 0.9 ±0.3, is compara- 
ble to values obtained for entanglement in the continuous 
variable regime [32[ and for entanglement of the discrete 
internal states of trapped ions [3, [4( . 
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FIG. 3: Concurrence C as a function of the storage time 
r. The solid line is obtained from Eq. J5]) assuming the 
fitted exponential decays, given in the inset, of the individual 
parameters p c and g±2 measured independently. The dotted 
line corresponds to Co. 



retrieval efficiency over the past year [l6|, [TtJ • Finally, by 
way of the conditional readout efficiency r] = 45 ± 10% 
for mapping of quantum states of the U, D ensembles to 
the fields 2^,2^, we estimate the density matrix pu,D 
and the concurrence Cjj,d = 0.9 ± 0.3 for the collective 
atomic state. We emphasize that Cu,d is an estimate 
determined from the model developed in Ref. [14| where 
the fields at the output of the MOT consist of a two-mode 
squeezed state plus background fields in coherent states. 

Turning then to a characterization of the decay of en- 
tanglement with storage time r, we present in Fig. [3] 
measurements of concurrence C(r) for fixed excitation 
probability p = 1.6 x 10 -3 corresponding to g±2 = 30 at 
r = 200 ns. C > for r < 20 /is, providing a lower 
bound for the lifetime of entanglement of the ensembles 
corresponding to 4km propagation delay in a fiber. 

To investigate the dynamics in Fig. 02 the inset shows 
the decay of the average gi2 and conditional probabil- 
ity p c for the ensembles taken independently. Such local 
decoherence has been investigated as the result of inho- 
mogeneous broadening of the Zeeman ground states due 
to residual magnetic fields [H, Ell, [31]. Our current mea- 
surement shows the effect of this local decoherence on 
the entanglement of the joint system of the ensembles. 
For this purpose, our measurements of C are superposed 
with a line C(r) given by Eq. ([5]), where the fitted expo- 
nential decay for p c (r), #12 (V) (with similar decay ~ 13 
[is) and the overlap £ determined in Fig. [2] are employed. 
The agreement evidenced in Fig. [3] confirms the principal 
role of local dephasing in the entanglement decay. 

In conclusion, we have reported a detailed study of the 
behavior of entanglement between collective excitations 
stored in atomic ensembles, including the dependence of 
the concurrence on the degree of excitation and the quan- 
titative relationship of local decoherence to entanglement 



[1] P. Zoller et a/., Eur. Phys. J. D 36, 203 (2005). 
[2] H.J. Briegel et al, Phys. Rev. Lett. 81, 5932 (1998). 
[3] H. Haffner et a/., Appl. Phys. B 81, 151 (2005). 
[4] C. Langer et al, Phys. Rev. Lett. 95, 060502 (2005). 
[5] B. Julsgaard, A. Kozhekin, E.S. Polzik, Nature 413, 400 
(2001). 

[6] L.M. Duan et al, Nature 414, 413 (2001). 

[7] C.W. Chou et a/., Nature 438, 828 (2005). 

[8] R.H. Dicke, Am. J. Phys. 49, 925 (1981). 

[9] C. Cabrillo et a/., Phys. Rev. A 59, 1025 (1999). 
[10] S. Bose et a/., Phys. Rev. A 83, 5158 (1999). 
[11] C.W. Chou et al, Science 316, 1316 (2007); published 

online 5 April 2007 (10.1126/science.ll40300). 
[12] J. Laurat et a/., New J. Phys. 9, 207 (2007). 
[13] Y.-A. Chen et a/., arxiv:0705:1256. 
[14] C.W. Chou et a/., Phys. Rev. Lett. 92, 213601 (2004). 
[15] V. Balic et a/., Phys. Rev. Lett. 94, 183601 (2005). 
[16] J. Laurat et a/., Opt. Express 14, 6912 (2006). 
[17] J.K. Thompson et al, Science 313, 74 (2006). 
[18] D. Felinto et a/., Phys. Rev. A 72, 053809 (2005). 
[19] H. de Riedmatten et al, Phys. Rev. Lett. 97, 113603 
(2006). 

[20] L. Jiang, J.M. Taylor, M.D. Lukin, Phys. Rev. A 76, 
012301 (2007). 

[21] Z.-B. Chen et al, quant-ph/0609151 

[22] B. Zhao et al, Phys. Rev. Lett 98, 240502 (2007). 

[23] W.K. Wooters, Phys. Rev. Lett. 80, 2245(1998). 

[24] D.N. Matsukevich, A. Kuzmich, Science 306, 663 (2004). 

[25] C.W. Chou, Towards a quantum network with atomic en- 
sembles, PhD Thesis, California Institute of Technology. 

[26] S.J. van Enk, N. Lutkenhaus, H.J. Kimble, Phys. Rev. A 
75, 052318 (2007). 

[27] C.H. Bennett et al, Phys. Rev. A 54, 3824 (1996). 

[28] D. Felinto et a/., Nature Physics 2, 844 (2006). 

[29] An alternative approach is to determine the suppres- 
sion of two-photon events relative to the square of the 
probability of single photon events for the fields 2, h = 
P11/ (pioPoi)- h < 1 is a necessary condition for entangle- 
ment 0[ll|,[l2[. Here, for g 12 = 60, h = 0.06. 

[30] J. Eisert, C. Simon, M. Plenio, J. Phys. A: Math. Gen. 
35, 3911 (2002). 

[31] D. Matsukevich et a/., Phys. Rev. Lett. 96, 033601 
(2005). 

[32] J. Laurat et al, Phys. Rev. A 71, 022313 (2005). 



